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ABSTRACT 


The  degree  to  which  the  ballistic  coefficient  ( 0  )  can  be  estimated 
from  reentry  metric  data  can  be  meaningfully  analyzed  with  a  sophisti¬ 
cated  computer  program  which  realistically  models  the  problem  and 
performs  a  proper  error  analysis  of  the  estimation  procedure .  Such 
a  program  has  been  developed  (at  AFETR)  and  is  being  used  to  isolate 
the  primary  error  sources  in  the  0  estimation  task.  This  computer 
program  is  additionally  used  to  determine  what  success  in  0  estima¬ 
tion  can  be  achieved  with  foreseeable  instrumentation  accuracies. 

Results  are  included  that  illustrate  the  effects  on  0  estimation 
from  the  following  variations:  1)  instrumentation  measurement  type  and 
accuracy;  2)  relative  geometry  between  trajectory  and  observer; 

3)  arc  length  and  minimum  altitude  of  tracking;  and  4)  the  magnitude 
Of  0  . 

The  central  thesis  of  the  analysis  presented  in  this  report  is  that 
dynamic  constraints  can  be  used  to  advantage  in  the  problem  of  estima¬ 
ting  the  trajectory  and  ballistic  coefficient.  The  least  squares,  con¬ 
strained  solution  takes  advantage  of  the  exercise  of  tracking  geometry 
over  a  sufficiently  long  arc  to  estimate  the  trajectory  parameters  and 
0  .  For  a  given  set  of  tracking  instrumentation,  the  enforcement  of 
dynamic  constraints  throughout  the  entire  trajectory  should  provide 
the  maximum  available  information  to  the  estimation  process. 
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I.  INTRODUCTION 


One  of  the  parameters  of  interest  in  current  post  flight 
reduction  of  radar  data  obtained  during  a  reentry  event  is  the  ballis¬ 
tic  coefficient  history  of  all  objects  tracked  by  the  observing  instru¬ 
mentation.  The  ability  to  successfully  determine  the  ballistic 
coefficient  0  from  the  measurement  data  is  influenced  by  the  data 
processing  techniques  employed.  Conversely,  a  given  accuracy  in 
the  0  estimate  can  be  achieved  with  minimal  instrumentation 
accuracy  if  the  most  efficient  data  processing  is  used  in  reducing 
the  measurement  data. 

This  paper  presents  the  results  of  an  error  analysis  of 
0  determination  capabilities  which  appraises  what  can  be  realistically 
achieved  with  foreseeable  coherent  radar  instrumentation,  and  what 
sources  of  error  can  significantly  degrade  such  estimates.  This 
study  was  conducted  using  the  recently  developed  BEAR*(Beta  Error 
Analysis  Routine)  program,  which  generates  a  trajectory  and  measure¬ 
ment  set  and  performs  a  minimum  variance  error  analysis  of  estimated 
variables,  which  in  this  case  are  the  trajectory  position  and  velocity  and 
the  ballistic  coefficient.  This  program  enforces  dynamic  constraints  in 
the  minimum  variance  error  analysis.  The  central  idea  of  this  analysis 
is  that,  for  given  instrumentation  accuracy  and  a  specified  tracking 
geometry,  the  use  of  the  dynamic  constraints  should  provide  maximum 
information  to  the  0  estimation. 

Interest  in  ballistic  coefficient  estimates  is  divided  into  two 
altitude  regimes;  the  free  molecular  flow  region  down  to  the  transition 
altitude,  and  the  lower  altitude  regime  after  continuum  flow  has  been 
established.  Therefore  the  analysis  which  follows  is  divided  into  an 
analysis  of  what  can  be  achieved  in  each  of  these  regions.  In  each  of 
these  flow  regimes  the  ballistic  coefficient  can  be  regarded  as  constant. 

*  BEAR  is  a  modified  version  of  the  ORAN  [  1  ]  program. 
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This  study  was  initiated  to  examine  the  degree  to  which 
estimation  of  the  ballistic  coefficient  can  be  improved  by  adding  a 
coherent  signal  processing  capability  to  the  advanced  Range  Instru¬ 
mentation  Ships.  This  addition  would  enable  the  ship  to  obtain 
independent  doppler  data  together  with  range,  azimuth  and  elevation 
measurements.  The  experiments  for  which  results  and  analysis  are 
presented  in  subsequent  sections  of  this  paper  were  designed  assuming 
the  sensor  to  be  a  ship  at  sea. 

n.  DESCRIPTION  OF  THE  ERROR  ANALYSIS  PROGRAM 

The  BEAR  program  is  an  IBM  360  program  for  computing 
the  effects  of  random  and  systematic  errors  on  minimum  variance 
orbit  determinations.  Systematic  errors  in  the  instrumentation  survey 
and  dynamic  model  may  be  considered  in  the  form  of  either  adjusted  or 
unadjusted  parameters,  with  the  effects  of  the  latter  broken  down  into 
individual  error  sources.  The  program  computes  the  effects  of  the 
unmodeled  parameters  on  both  the  orbit  and  recovered  parameters, 
with  the  orbital  effects  propagated  from  epoch  to  any  desired  prediction 
time. 

Geopotential  parameters  plus  one  drag  parameter  may  be 
considered.  Partial  derivatives  with  respect  to  these  parameters 
are  computed  by  numerical  integration. 

The  ephemeris  of  position  and  velocity  is  computed  analytically 
by  an  orbit  generator  subroutine  [2]  .  In  this  subroutine  the  dynamic 
partial  derivatives  relating  position  and  velocity  at  any  time  to  orbit 
parameters  at  epoch  are  also  computed  analytically,  as  are  the  first 
order  secular  perturbations  due  to  Jg  and  atmospheric  drag.  Some 
of  the  mathematics  [3]  will  be  presented  here,  in  a  compact  form,  to 
describe  what  is  being  computed. 
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First,  hypothesize  the  linearized  equations  in  vector  and 
matrix  notation  as 


y=By  +  U  6  +  e  ^ 

nxl  nxppxl  nxu uxl  nxl 

where 

y  -  observation  vector  from  linearization;  i.  e. ,  the 
difference  between  the  actual  observations  and 
the  approximate  value 

y  -  adjusted  parameter  vector  corrections  (consists 
of  corrections  to  the  orbital  elements,  survey, 
instrumentation  biases,  etc. ) 

B  -  matrix  of  partial  derivatives  of  the  observations  (O) 
and  the  adjusted  parameters 

6  -  unadjusted  or  unmodeled  parameter  vector  errors 
(could  consist  of  geopotential  uncertainties,  survey 
errors,  etc.) 

U  -  matrix  of  partial  derivatives  of  the  observations  and 
the  unadjusted  parameters 

e  -  vector  of  true  observation  errors.  Assume  Var  (e)  =  Z 
a  diagonal  matrix.  n* 


The  standard  least  squares  solution  for  the  adjusted  parameters, 
when  6=0,  is 


Y  =  (B 


y 


(2) 


where 


Var  (y)  =  (BT 


(3) 


Equation  (3)  gives  the  standard  result  for  the  accuracy  of  the  actfusted 
parameters.  Conventionally,  the  square  roots  of  the  diagonal  elements 
of  (B  Z  B)“  are  the  standard  deviations  or  sigma  values  for  the 
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adjusted  parameters.  Now  when  5  4  0  ,  it  can  be  shown  that 

Var  (y)  =  (BTL"1B)  1  +  DWuDT  (4) 

where 

D  =  (BTS_1B)  1BTZ"1U=  —  (5) 

as 

Wu  is  the  variance-covariance  matrix  of  the  unadjusted  parameters  6  , 
and  is  taken  as  a  diagonal  matrix.  The  square  roots  of  the  diagonal 
elements  of  (4)  now  include  a  contribution  due  to  the  unadjusted  effects. 

The  dynamic  constraints  are  contained  in  the  B  matrix, 
which  is  computed  as  the  matrix  product 


0Oi  3Xi 

"  9"x7  TT0 

and  o  denote  the  ith  time  point  and  epoch 
respectively.  The  dynamic  constraints  are  contained  in  the  second 
term  of  equation  (6)  . 

For  this  study,  the  BEAR  computer  program  was  used  to 
provide  a  minimum  variance  estimation  of  the  ballistic  coefficient 
simultaneously  with  the  trajectory,  where  dynamic  constraints  were 
enforced.  Computer  runs  were  made  to  examine  the  effects  on  0 
estimation  from  the  following  variations: 


B  = 


30. 
3r  o 


where  the  subscripts  i 


(6) 


(1)  Instrumentation  measurement  type  and  accuracy 

(2)  Relative  geometry  between  trajectory  and  observer 

(3)  Trajectory  arc  length  and  minimum  altitude  of  tracking 

(4)  The  magnitude  of  the  ballistic  coefficient 

(5)  Instrumentation  measurement  biases 

(6)  Survey  uncertainties 

(7)  Uncertainties  in  the  location  of  the  dynamic  center  of 
mass  of  the  earth 

(8)  The  addition  of  a  second  sensor. 

The  results  of  these  studies  are  discussed  in  the  next  section. 
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III.  ERROR  ANALYSIS 


In  the  high  altitude  free  molecular  flow  region,  the  atmos¬ 
pheric  drag  force 


F n  =  -Ip  mg  ,  with 

D  P 


8  -  W  lb 

CnA  ft2  ’ 


(7) 


is  quite  small  [4  ] .  Attempts  to  measure  these  minute  decelerations  at 
each  point  in  time  and  thereby  determine  /3  on  a  point-to-point  basis 
translate  into  unattainable  instrumentation  accuracy.  However,  a  long 
arc  solution  which  exploits  dynamic  constraints  and  variable  geometry 
can  see  a  cumulative  effect  which  is  not  detectable  with  point-to-point 
observations. 

As  an  object  makes  the  transition  from  free  molecular  to 
continuum  flow,  the  ballistic  coefficient  increases.  For  purposes  of 
discussion,  it  is  reasonable  to  assume  that  an  object  which  has  hyper¬ 
sonic  continuum  flow  /3  of  500-1000  will  have  a  free  molecular  flow 
ballistic  coefficient  of  about  150-200.  In  all  cases  presented  here  no 
a  priori  information  was  assumed  for  the  trajectory,  which  is  estimated 
simultaneously  with  the  ballistic  coefficient. 

Figure  1  illustrates  the  results  of  a  BEAR  analysis  with 
the  parameter  a  ft  chosen  as  the  figure  of  merit  in  the  following 
curves.  The  true  ballistic  coefficient  for  this  curve  is  160  and  the 
reentry  parameter  at  400K  feet  are  V  g  =  34°  and  =  2x10^  ft/sec. 
A  single  doppler  radar  with  the  indicated  measurement  capability, 
located  approximately  50  miles  downrange  of  impact,  is  the  observing 
instrument.  These  accuracies  were  chosen  to  represent  reasonable 
state-of-the-art  capability  in  the  future.  Moreover  ,  very  large 
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FIGURE  1 

High  Altitude  p  Estimation 
Using  Dynamic  Constraints 


values  were  assumed  for  the  a  priori  cr  ^  to  simulate  a  lack  of 
knowledge  of  the  value  of  0.  By  this  means  the  estimate  is  based 
entirely  on  measurement  information  collected  during  the  event.  This 
curve  has  been  extended  to  lower  altitudes  because  the  transition  point 
is  dependent  on  body  shape  according  to  the  relation  .  01  so 

that  the  lower  cutoff  point  is  variable,  where  X  is  an  atmospheric 
mean  free  path  and  D  is  the  body  diameter.  Initial  acquisition  of  the 
object  occurs  approximately  6.  5  minutes  before  impact  at  a  range  of 
103  miles.  The  measurements  are  summarized  in  Figure  2. 

Under  these  conditions,  the  curves  show  that  it  is  unlikely 
that  0  can  be  resolved  to  within  a  factor  of  two  for  objects  of  this  class 
( 0  cont  ~500  "  1000 )  before  the  transition  flow  regime  occurs.  If  another 
identical  radar,  located  some  50  miles  uprange  and  displaced  from  the 
trajectory  is  available  to  assist  in  the  discrimination,  Figure  1  shows 
that  the  situation  is  only  slightly  improved  in  the  high  altitude  region. 

The  ability  to  estimate  0  is  not  directly  a  function  of  the 
ability  to  measure  position  and  velocity.  Table  1  shows  the  position  and 
velocity  standard  deviations  at  points  of  interest  on  the  trajectory  corres¬ 
ponding  to  Figure  1.  Notice  that  even  though  the  trajectory  is  improved 
by  a  factor  of  two,  the  0  estimate  improves  by  approximately  20%.  The 
dynamic  constraint  filter  is  greatly  enhanced  by  the  doppler  measurement 
of  one  component  of  velocity,  so  the  optimum  configuration  is  a  sensor 
located  to  maximize  the  doppler  component  of  velocity  over  the  0  estima¬ 
tion  interval.  Table  2  gives  a  comparison  of  the  ability  to  estimate  0  and 
the  trajectory  with  and  without  an  R  measurement,  for  two  values  of  0  . 
Even  though  the  doppler  measurement  does  not  greatly  affect  the  trajectory 
determination,  it  directly  influences  the  0  estimate.  In  fact,  studies 
have  shown  that  for  a  given  geometry  and  instrumentation  accuracy, 
very  nearly  scales  with  . 

A  determination  of  the  ballistic  coefficient  at  lower  altitudes 
in  continuum  flow  regime  is  of  interest  primarily  for  intelligence  purposes 
rather  than  discrimination  information  on  uncooperative  objects. 


FIGURE  2 


Measurement  History  For 
P  =  160  Trajectory. 


Time  From  Epoch  (minutes) 

1  2  3  4  5  6 

1 - 1 - 1 - 1 - 1 - r 

Altitude  (feet  x  10®  ) 


TABLE  1 


Single  Ship 

Two  Ships 

Time  From  Epoch  (min) 

°p  <ft> 

<JV(-£L-) 
v  sec 

ap 

ffv 

0 

8.67 

.03 

4.66 

.01 

1.0 

7.25 

.03 

3.88 

.01 

2.0 

5.71 

.03 

3.02 

.01 

3.0 

4.29 

.03 

2.24 

.01 

4.0 

2.99 

.03 

1.52 

.01 

5.0 

2.05 

.03 

1.02 

.01 

6.0  (320  K  ft) 

2.03 

.03 

1.09 

.01 

6.2  (197  K  ft) 

2.16 

.03 

1.18 

.01 

6. 35  (cut  off  measurements) 

2.27 

.08 

1.24 

.06 

3.7  K  ft. 


TABLE  2 


/ (3  a  priori  =  2. 0 

aR  =  -1 

• 

No  R  Data 

0 

a  a 

v  p 

ye 

a 

V 

a 

P 

160 

.12 

.03  ft  2.0  (ft) 

1.9 

.06 

3.1 

sec 

1600 

.23 

.02  1.5 

1.97 

.06 

3J 

The  ability  to  estimate  0  below  300K  feet  as  a  function  of 
acquisition  altitude  is  shown  in  Figure  3  for  two  different  values  of  the 
ballistic  coefficient.  The  power  of  dynamic  constraints  is  illustrated 
in  Table  3,  which  compares  a  ^  /  0  as  determined  by  the  BEAR 
program  with  values  obtained  by  ten  second  filtering. 

The  ability  to  determine  0  is  a  function  of  instrumentation 
location  relative  to  the  trajectory  plane.  Table  4  compares  a ^  /0 
determinations  for  radar  locations  50  miles  downrange,  crossrange, 
and  at  a  45°  angle  with  respect  to  the  trajectory  plane. 

As  a  final  topic,  an  example  of  the  effects  of  unadjusted 
parameters  on  0  determination  will  be  presented.  A  single  station 
solution  is  considered  with  the  standard  deviation  of  the  errors  shown 
in  Table  5. 

If  propagated  as  uncorrelated  errors,  they  contribute  to  the 
error  in  0  as  shown  in  Table  6.  If  these  error  sources  exist  in  the 
instrumentation  it  would  be  necessary  to  adjust  for  these  parameters 
along  with  0  .  Table  6  shows  the  results  of  an  attempt  to  adjust  the 
errors  and  the  residual  effect  on  a ^  . 

Survey  errors  for  a  single  station  need  be  included  only  if 
it  is  necessary  to  express  the  measured  trajectory  relative  to  a  known 
coordinate  system.  For  example,  this  will  be  the  case  when  it  is  desired 
to  determine  the  launch  point  for  intelligence  purposes.  For  0  deter¬ 
mination  it  is  necessary  to  estimate  the  trajectory  relative  to  a  single 
station,  hence  only  the  error  in  the  dynamical  center  of  the  earth 
relative  to  the  ship  need  be  considered.  An  analysis  similar  to  that 
of  Table  6,  where  the  relative  uncertainty  in  the  radial  distance  to  the 
dynamical  center  was  assumed  to  be  150  feet,  has  shown  that  this  effect 
is  less  than  one  percent  on  /  0  .  Survey  errors  could  also  be  signi¬ 
ficant  when  it  is  desired  to  combine  the  measured  trajectory  with  data 
from  another  station.  Consider  again  the  two-station  solution  of 
Figure  1,  where  the  second  radar  improves  the  0  determination  by 
approximately  twenty  percent.  If  an  error  of  2500  feet  is  assumed 
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TABLE  3 


a  priori  c/^  y/3  =10 


V*3 

v*3 

Dynamic 

10  Sec 

Constraints 

Filter 

150  K  ft, 

0  =  800 

.  55  (Fig.  3) 

1.84 

250  K  ft, 

/3  =  800 

.42  (Fig.  3) 

10 

300  K  ft, 

/3  =  160 

1.5  (Fig.  1) 

10 

1 

TABLE  4 

p  =  800  lb/ft2 

% 

Downrange 

.20 

1.75 

.02 

45°  Angle 

.32 

1.2 

.02 

Crossrange 

.35 

1.35 

.02 

This  table  again  illustrates  the  power  of  the  R  measure 
ment  in  determining  (i  . 
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Acquisition  Altitude  (feet  x  1(T  ) 


•It 


.iTi 


1 


TABLE  5 


a  priori  o 

Adjusted  a 

1. 

R  Bias 

50  ft 

21.3  ft 

2. 

A  Bias 

.  00066  mil 

.  00004  mil 

3. 

E  Bias 

.  00066  mil 

.  00011  mil 

4. 

E  Refraction 

.  1  mil 

.05  mil 

5. 

X  Survey 

2500  ft 

2144.0  ft 

6. 

Y  Survey 

2500  ft 

1175.0  ft 

TABLE 

ap/P 

No  Biases 

Ballistic 

Coefficient  .  208 

_6 

Unadjusted  Adjusted 

Biases  Biases 

18.95  .431 

Individual  Effects  of  Unmodeled  Parameters  on  a ^ 

Unadjusted 

Contributions 

Parameters 

To  cr^  //3 

R  Bias 

.75 

A  Bias 

.37 

E  Bias 

18.6 

E  Refraction 

.4 

X  Survey 

2.9 

Y  Survey 

1.  5 
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along  two  aces  for  the  uncertainty  in  the  relative  location  of  one  of  the 
ships,  these  bases  will  completely  negate  the  improvement  in  /3  esti¬ 
mation  gained  by  the  presence  of  the  second  ship. 

It  is  of  particular  interest  to  the  Aris  ship  analysis  to 
determine  what  precision  in  R  is  required  to  improve  the  /3  determi¬ 
nation  capability.  Table  7  presents  a  matrix  of  <J  p  /  P  ,  estimation 
capability  as  a  function  of  R  precision  and  ship  location  relative  to 
impact.  The  present  capability  with  no  R  is  also  included  for  compari¬ 
son.  This  data  shows  that  little  benefit  is  derived  from  R  measurements 
with  c >  .2  ft/sec  ,  where  the  ship  is  located  within  a  45°  cone  with 
respect  to  the  trajectory  plane. 


TABLE  7 

/3  =  800  lb/ft2 

CT/3 

//3  a  priori  =  1.0 

aR 

2.0 

Downrange 

45°  Angle 

Crossrange 

1.0 

1.0 

1.0 

.2 

.65 

.65 

.92 

.  1 

.20 

.32 

.35 

No  R 

1.0 

1.0 

1.0 
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IV.  CONCLUSIONS 


In  both  free  molecular  and  continuum  flow  the  use  of  dynamic 
constraints  appears  to  have  potential  for  0  estimation. 

For  slender  bodies  where  the  atmospheric  drag  force  is 
relatively  small,  and  in  the  high  altitude  regions  of  the 
atmosphere,  the  dynamic  constraint  technique  appears  to 
offer  significant  advantage. 

For  the  cases  examined,  range,  azimuth,  and  elevation 
primarily  determine  the  trajectory,  while  the  inclusion 
of  an  R  measurement  greatly  increases  the  capability 
to  estimate  0  . 

Because  of  the  importance  of  R  in  the  determination  of  0  , 
it  is  desirable  to  align  the  sensor  with  the  trajectory  plane. 

Bias  errors  in  the  measurements,  especially  elevation 
angle  errors  in  this  case,  seriously  degrade  0  estimation 
if  they  cannot  be  removed  by  calibration. 

Data  reduction  techniques  that  can  simultaneously  adjust  for 
the  trajectory,  0  and  instrumentation  biases  appear  to 
offer  the  possibility  of  significant  improvement  in  0 
estimation. 

Survey  errors  for  a  single  sensor  are  significant  only  when 
transforming  the  measurements  from  the  sensor  referenced 
coordinate  system;  otherwise  they  do  not  influence  the  0 
determination. 

Uncertainties  in  the  location  of  the  dynamical  center  of  the 
earth  has  negligible  effect  on  0  estimates  obtained  using 
dynamic  constraints. 

The  addition  of  a  second  sensor  near  the  impact  point  contri¬ 
butes  on  the  order  of  twenty  percent  to  the  0  estimation 
capability.  However,  relative  survey  errors  can  negate  this 
improvement. 
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